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Abstract: This paper provides a review of some recent is- 
sues on the Mixmaster dynamics concerning the features of its 
stochasticity. After a description of the geometrical structure 
characterizing the homogeneous cosmological models in the 
Bianchi classification and the Belinsky-Khalatnikov-Lifshitz 
piecewise representation of the types VIII and IX oscillatory 
regime, we face the question regarding the time covariance of 
the resulting chaos as viewed in terms of continuous Misner- 
Chitre like variables. Finally we show how in the statistical 
mechanics framework the Mixmaster chaos raises as semiclas- 
sical limit of the quantum dynamics in the Planckian era. 

PACS number(s): 04.20. Jb, 98.80. Dr 



I. INTRODUCTION 

As well-known the Standard Cosmological Model 
(SCM) finds its theoretical basis in terms of a homo- 
geneous and isotropic universe obtained as high sym- 
metry solution of the Einstein's equations, the so called 
Friedmann- Robertson- Walker model (FRW). Such repre- 
sentation of our actual universe possesses a clear degree 
of reliability due to its good general agreement with re- 
spect to the observed phenomenology (in particular the 
strong isotropy of the Cosmic Background Radiation as 
well as the consistency between the predictions on the 
primordial nucleosynthesis of the light elements and the 
experimentally observed abundances), nevertheless there 
are some important general aspects to be taken into ac- 
count. 

In first place none theoretical principle led us to ex- 
clude that in the very early phases of its evolution the 
universe had been characterized by a higher degree of in- 
homogeneity and anisotropy and only in a later stage un- 
derwent an isotropization process as natural consequence 
of its dynamics and/or by the action of some physical 
mechanism in classical as well as quantum regime (to 
which the causality notion plays a crucial role); indeed 
the instability of the FRW solution toward the cosmo- 
logical singularity implies that a more general behavior 



should characterize the very early evolution as soon as 
we make allowance for small perturbations responsible 
for the actual clumpyness. 

Furthermore the SCM in its original formulation con- 
tains non-trivial internal inconsistencies which require 
to be explained through appropriate modifications of 
the underlying theory. Among this undesired theoreti- 
cal facts take particular importance the so called "hori- 
zon paradox" and "flatness problem" to which should be 
added the absence of an appropriate model for the large- 
scale structures formation able to reproduce t he observ ed 
distribution of matter in the actual universe 
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All these considerations make clear the deep interest 
in studying more general classes of solutions of the Ein- 
stein's equations in order to individualize dynamical be- 
haviors which could constitute a more suitable frame- 
work than the simple FRW model for the construction 
of a completely self-consistent cosmology, with particu- 
lar reference to the universe evolution in proximity of the 
initial "Big Bang" |4|j3|j5|. 

In this paper (Sections I and II) we will discuss one of 
the most studied cosmological dynamical models in view 
of its possible implications in the history of our universe, 
the so called Mixmaster model It corresponds 

to the asymptotic evolution toward the cosmological sin- 
gularity of the type models VIII and IX in the famous 
Bianchi classification |l^. Another important theoret- 
ical property characterizing such model is that it is an 
appropriate prototype of the behavior of the general cos- 
mological solution of the Einstein's equations in the same 
asymptotic region to the initial singularity. Indeed the di- 
rect extension to the inhomogeneous case of the ideas and 
the formalism required by the treatment of the homoge- 
neous one, extension based on the dynamical decoupling 
characterizing near the singularity the different points of 
the space, leads to derive in natural way the asymptotic 
evolution to the "Big Bang" of a generic inhomogeneous 
cosmological model which constitutes surely one of the 
most important results until now obtained in Relativis- 
tic Cosmology [0,|9|. 

The description of the Mixmaster Model can take place 
in terms of two different, but indeed completely equiv- 
alent approaches: one is based on the direct analysis of 
the Einstein's equations, as in the original analysis due to 
Belinski, Khalatnikov and Lifshitz (BKL) (Section III), 
while the other corresponds to a Hamiltonian approach 
to the dynamics (Section IV) which had been introduced 
first by Misner M. 
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However not all anisotropic dynamics are compatible 
with a satisfactory SCM but, as shown in the early Sev- 
enties, under suitable conditions some can be represented 
as a FRW model plus a gravitational waves packet ( , 

Among the Bianchi classification, the types Vllljl] and IX 
appear as the most general ones: the former's geometry 
is invariant under the S0{3) group, like the closed FRW 
universe and its dynamics allows the line element to be 
decomposed as 



ds'^ = dso - Sfa 



(a)(b)^ifc 



(1) 



where dso denotes the line element of an isotropic uni- 
verse having positive constant curvature, Gj-^'"'^'' is a set 
of spatial tensors^ and S(^a){b){t) are amplitude functions, 
resulting small sufficiently far from the singularity. 
Since Belinski, Khalatnikov and Lifshitz Q derived the 
oscillatory regime characterizing the evolution near a 
physical singularity, a wide literature faced over the years 
this subject in order to understand the resulting chaos. 

The research developed overall in two different, but re- 
lated, directions: firstly, the dynamical analysis removed 
the limits of the BKL approach due to its discrete nature, 
providing satisfactory representations in terms of contin- 
uous variables (i.e. construction of an invariant measure 
for the system ||2^], [Q); secondly, a better characteri- 
zation of the Mixmaster chaos (in view of its properties 
of covariance) found non-trivial difficulties to apply the 
standard chaos indicators to relativistic systems. 
The various approaches prevented, up to now, to say a 
definitive word about the covariance of Mixmaster chaos, 
though important indications in favor of such covariance 
arise from Q, |Q, Q and js^js^, as above discussed 
in detail. 

A wide interest in understanding this intrinsic nature 
of the chaotic Mixmaster dynamics and its very early ap- 
pearance in the Universe evolution, to which we dedicate 
our analysis, lead to believe in the existence of a relation 
with the quantum behavior the system performs during 
the Planckian era. In fact the last aim of this paper 
(Section V) is to give a precise meaning to such relation 
by constructing the semiclassical limit of a Schoedinger 
approach to the canonical quantization of the Arnowitt- 
Deser-Misner (ADM) dynamics. 



^ All the considerations we will develop for the type IX apply 
also to the VIII one since, close to the singularity, they have 
the same morphology. 

^ These tensors satisfy the equations 



{a){6) 
ik ' 



(a)(b)k 



= 0, G' 



(a)(6)i 



= 0, 



in which the Laplacian is referred to the geometry of the 
sphere of unit radius. 



II. GEOMETRICAL STRUCTURE OF THE 
BIANCHI MODELS 

A. Group Transformations 

A (pseudo) Riemannian space, given a metric g and a 
differential structure (M, g), is homogeneous if it is orbit 
for an isometry group as a group invariance for the met- 
ric g. The groups of movements for the metric are said 
to be isometry groups and for homogeneous spaces they 
are transitive, in the sense that each point is equivalent 
under the action of the group and they are spatial sec- 
tions of the space-time. 

We are interested in the Lie algebras for Killing vector 
fields as generators of these groups of movements, in- 
tended as infinitesimal transformations generating rota- 
tions and translations. 
Let's consider a group of transformations 



x'^ -> = ft" (x, a) 



(2) 



over a space M, where {a°'}^^i ^ are r independent 
variables parameterizing the group and ao corresponds 
to the identical transformation 



(3) 



Let's then consider also the infinitesimal transformation 
corresponding to oq -t- <5a, close to the identity. 



= /A" (x.^ ao + Sa) w 



(4) 



say 



« x'^ + (x) Sa" = (1 + Sa^^a) , 



(5) 



where the r first order differential operators {^a} are de- 
fined as = in correspondence with the r vecto- 
rial fields with components {^^}. These are the Killing 
generating vector fields. 

The coordinate transformation can be definitely written 
as 



x'^ « (1 + Sa'^^a) 
and in a finite form 



x'^ 



(6) 
(7) 



where {0°} are r new parameters for the group. 
The vectorial generator fields form a Lie algebra, say a 
real r-dimensional vectorial space with basis {Ca}i closed 
with respect to commutation, in order to have the repre- 
sentation 



(8) 
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where C"^^^ are the structure constants for the Lie algebra. 
It is then natural to extend this formalism to the choice 
of {ca} as basis for the Lie algebra g over a group G with 

[ea,et]=C\,e,; (9) 

then let's define the symmetric quantity 

Jab = = Iba (10) 

and an internal product over the Lie algebra 

7afc = • = 7 (e„ e^) , 7 (X^^e^, Y^'e^) = ^labX'^Y'' . 

(11) 

The orbit of x is 

/g {x) - {/a (a:) I a e G} (12) 

as the set of all points that can be achieved from x under 
the group of transformations. The isotropy group in x is 

G^^{a£G\fa {x) = x} , (13) 

i.e. the subgroup of G which leaves x fixed. If xi,X2 
are on the same orbit, then Gxi and Gx2 are conjugate 
subgroups of G and then isomorphic. If now Gx = {oq} 
and fa [x) ~ M, G is diffcomorphic to M and the two 
spaces are identified. If g is a metric over M invariant 
under G, then it is definitively given by the inner product 
of the invariant vectorial fields of basis ea- 

The groups of non-Abelian transformations, as defined 
by (^), represent geometrically homogeneous spaces in 
three dimensions and give the section spatially homoge- 
neous of the spatially homogeneous space-times. 
Given a basis {ca} of the Lie algebra for the tridimen- 
sional Lie group, with structure constants G'^^j,, at any 
time the spatial metric is given by the spatially constant 
inner products 

Ra- e.b= gab{t) a, 6=1,..., 3 (14) 

which are six functions of the time variable only. This 
permits to rewrite the Einstein equations as ordinary dif- 
ferential equations, eventually associated to constraints 
functions of t, necessary to describe the matter behavior 
of the universe. 

In four dimensions one obtains homogeneous space-times: 
given a set of structure constants G"^^ over a basis {bq} 
of the four dimensional Lie algebra, the ten constants 

5a/3 = 6^6/5 a,/3=l,...,4 (15) 

determine univocally the signature of the Lorentz metric. 
As a consequence, Einstein's equations reduce to a sys- 
tem of algebraic equations for gap and G"^^, which could 
in principle not be solvable for each transformations' 
group, in view of the fact that for general formulations 
with matter fields there are more constants to take ac- 
count of. 



In any case one needs to consider only one representative 
group for each class of equivalence of the Lie groups. Re- 
ferring to three dimensions (only the spatial sector), the 
Bianchi classification determines definitely all symme- 
tries for tridimensional homogeneous spaces, analogously 
to the curvature (fc = —1, 0, +1) which distinguishes ho- 
mogeneous and isotropic spaces (FRW). 

B. Einstein's Equations in the Synchronous Gauge 

Let's us specify the general scheme outlined above in 
the specific approach of Luigi Bianchi in 1897 [|] and 
independently applied to cosmology by Belinski, Khalat- 
nikov and Lifshitz in 1969 0|. 

In the cosmological approach we consider a synchronous 
reference system, so that there is an unique time t syn- 
chronized all over the space. This can be obtained requir- 
ing the metric tensor gab to have the peculiar diagonal 
form such as 

500 = 1 go„ = 0, q; = 1,...,3 (16) 

in order to write the line clement 

ds^ ^ N^{t)dt'^ - df , (17) 

where N{t) is the lapse function (in this section we adopt 
the synchronous gauge = 1). Writing^ 

df^jaf3{t,x'>)dx°'dx'^, (18) 

t is the synchronous time and the tridimensional ten- 
sor "faf3 describes the metric spatial sector. The iden- 
tity goa = ensures the complete equivalence in the 
three spatial directions, while the unitarity of goo is 
always valid, provided an appropriate rescaling of the 
time coordinate. In this reference system the time lines 
{x'' = const.) are the geodesies of the four dimensional 
space-time 

dii^ 

— +Tliu''u^=n„=0, (z,fc,Z = 0,...,3) (19) 

being = ^ the tangent four vector to the world line 
x^ — const, with components it" = 1, it" = and or- 
thogonal to the hypersurfaces t — const.^. 



■^All over this section Latin indexes run over 0, . . . , 3 while 
Greek indexes 1, . . . , 3. 

''This geometrical construction is always available but not 
univocally defined: a coordinate transformation of the type 

( t' =t 

doesn't affect the time coordinate and defines the passage 
from a synchronous reference system to another one, fully 
equivalent. 
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Under this choice and the metric ( p7| ) , the mixed compo- 
nents of the Einstein's equations write 



inverse matrix is rj , being r]°''^r]cb = S^- Then follow the 
properties 



R' 



(22) 



Ri = -Pf - 77^1 (77'^^) = 8^G - ij^T ) (23) 



being 



dial: 

dt 



7 =1 lap I 



(24) 



Tap the momentum energy tensor for the system and Pap 
the tridimensional Ricci tensor obtained from the metric 

lapK 



In ( p3D the spatial and temporal derivatives are split, 
while from ( pT| ) one obtains the property for the met- 
ric determinant to be zero for a certain time (Landau- 
Raichoudhury theorem)^ nevertheless the singularity 
with respect to the temporal parameter is a physical one. 



C. Tetradic Representation of the Metric Tensor 

The choice of a tetradic basis of four linearly indepen- 
dent vector fields (depending from the system symme- 
tries) permits to project all quantities in a very useful 
way to obtain new simpler equations to be satisfied. 
Consider on each point of the space-time the basis of 
four linearly independent contravariant vectors e^^^ , (a = 

1, . . . ,4), where the bracketed index is the tetradic one, 
while i is tensorial, and the covariant set e(a) i = gike^a)^ 
being gik the metric tensor. Define also the reciprocal 
set e|"''eJ^^ = S'^ (orthogonality condition). It is easy to 



check the validity of 

The definition is complete with 

e(a)e(6) I = Vab : 



(25) 



(26) 



where rjab is a symmetric matrix with signature (H 

) (with the orthonormality of (|25|)); the corresponding 



V{a){b)e, 



(a) 



Hb)z 



'/ t^{a)i — 

(a) _ 



(27) 
(28) 
(29) 



For a generic vector (Aj) or tensorial (T^) field the 
tetradic components are, in general 



A' 



and, respectively. 



^ (a)(6) 

T 



(a)^\b)-^i3 - ^\a)Ti{b) , 
(a) (b)rri 



(30) 

(31) 
(32) 

(33) 
(34) 



The tetradic indexes can be raised or lowered with the 
use of the tensors ?7(a)(6) and rj^"^^''^ , while the contraction 
gives then a result independent of the indexes nature. 
By (26) and ( p7|) one obtains also 



gik = e(a)»efe"' = Vabefe}^' (35) 



so that the line element becomes 

ds' = Vab {ei'^dx^) {ei'Ux'^) . 



(36) 



The choice of rjab permits to split the tetradic basis in 
one temporal and three spatial vectors. Nevertheless the 
expressions dx^°'^ = ef^dx^ are not, in general, exact 
differentials of functions of the coordinates. 



D. Directional Derivative 

The contravariant set e(a) of tangent vectors leads to 
the natural covariant derivative definition 



S(a) 



-{a) Q^i 



(37) 



so that the derivative of a generic scalar field along the 
direction a is 



^The indexes of these two tensors are raised and lowered 
using such reduced metric. 

®This fictitious singularity disappears with a change of refer- 
ence system and is related to the intersection of the geodesies 
belonging to any arbitrary set with the envelope surfaces, 
say to the procedure used to build a synchronous reference 
system. 



The general extension of such definition is 

, d 



^{a),(b) — e)h\Tr^A 



9 i , 
e'(h\Tr-e-',.Ai 



in order to rewrite 



'(b) 



(38) 



(39) 
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Let's introduce the Ricci's rotation coefficients ^ahc 

labc = (^(a)i-k&\b)(i\c) (41) 

and their hnear combinations 

Aafcc = labc — lacb = {<^{ayi;k " ^(^a)k-i) ^\b)^\c) = 

= (e(Q)i,fc - e(Q)fc,i) e(fc)e(^^) (42) 
in which is has been used the identity 

dA, dAk 



R(a)(b)(c)(d) = Rniikie'(lflb)e^cf[a) 



Ai-k ~ Ak-i 



(43) 



Expression ([42|), in which the regular derivatives are sub- 
stituted by the covariant ones, is invertible 



labc — 2 {^abc + ^bca ^ ^cab) 



From the identity 

= V{a){b), 



e-{a)ke{b) 



one gets the symmetry properties 



labc 



-Ibac 



X, 



.be 



-X 



acb ■ 



(44) 



(45) 



(46) 



Now the formahsm is ready to find the values of the struc- 
ture constants which leave the metric invariant under the 
homogeneity constraint. 

The basis 6(3) permits to express the Lie parentheses as 



[e(a),e(b)] = C^^l^^^kfic) 



(47) 



where the coefficients C^^^jj-^^ are the 24 (in four dimen- 
sional space) structure constants for the group of trans- 
formations, antisymmetric with respect to the lower in- 
dexes; it is easy to obtain the explicit relation 



C 



ic) 
(a)(b) 



(c) _ (c) 

T (b)(a) 7 (a)(6) 



(48) 



E. Ricci and Bianchi Identities 

By the Riemann tensor -R™^; , the Bianchi identity be- 
comes 



Ai-k:i ~ Ai-i-h — AfnR"ii., 



kl 



(49) 



for a generic four vector Ai and, applied to the tetradic 
basis, 

e{a)i-kd ~ ^{a)i;l;k = e^)Rrmkl ■ (50) 

Projecting this expression on the basis itself one obtains 



-7(a)(6)(c), (d) +7(a)(b)(<i),(, 
"7(a)(6)(/) 



(f) (f) 
l(cf\d)-l{d)''\c) 



if) 

-lU)(a)(c) 1(b)'" id) 
( f) 

-7(/)(a)(d) l{b)'"{c) 



(51) 



In view of the re-expression of the interesting quantities 
in terms of the Ricci's rotation coefficients l[a)(b)(c) f^^d 
subsequently in terms of the structure constants, it is 
easy to see the great simplification in the formulas, pro- 
vided a space for which exists a set of such constants 
^'^(6)(c)' ^^^'^ analysis made by Bianchi for the homo- 
geneous spaces leads to [0 . 

The length element before the transformation is 



lap 



x^,x^\ dx°'dx'^ 



(52) 



that, under a general change of coordinates, transforms 
to 



dl^ 



7„/3 {x'\x'^,x'^)dx""dx' 



(53) 



where the functional form jap has to be the same un- 
der the homogeneity constraint. Such requirement, in a 
uniform non Euclidean space, leads to invariance of the 
three independent differential forms in ( ^6| ) which are not 
exact differential of any function of the coordinates and, 
in tetradic form, are 



e'^^dx" . 

The spatial invariant line element rewrites 

dp = ryafc [el^^dx'^) {4''dx^) 
so that the metric tensor becomes 
la/3 = 77a6ei°^e^''^ , 



(54) 



(55) 



(56) 



maintaining for rjab the definition (Pq), as a function of 
the time variable only. 

The symmetry properties for the space determine the 
specific choice of the basis vectors which, in general, are 
not orthogonal and consequently the metric r]ab is not 
diagonal. 

In pure specific case the relations between the three vec- 
tors are 



(57) 



6(1) = 


1 


e(2) 


Ae(3) 


V 






6(2) = 


1 






V 






6(3) = 


1 




Ae(2) 


V 







where v is the product 



=(a) 



,(1) 



(58) 
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and where it is natural to interpret eja) and e^") as the 

Cartesian vectors of components e"^^ and eL"'' respec- 
tively. 



Multiplying (65) by ej^^ gives the uniformity conditions 
over the space 



The metric tensor determinant (56) takes the value 

7 — rjv^ (59) 

being 77 the determinant of the matrix rjab- 

The space invariance as in ( p3[ ) is equivalent to the in- 

variance of (IsJ), then 



ifW ifW 



(66) 



The expression on the left hand side corresponds to the 
definition of A'^^^ (^), then constants. 
The antisymmetry property holds, see ( ^ ) or (^8|), with 
respect to the lower indexes 



e(f ) {x) dx" = el^^ {x') dx" 



(60) 



(67) 



where eL"'' on both sides of the equation are the same 
functions of the old and new coordinates respectively. 
Making some algebra one gets 



Such relations can be rewritten in a compact form in 
terms of the linear operators 



d 



(68) 



dx^ 



(61) 



which is a system of differential equations to determine 
the functions x'^ (x) on the given basis. The integrabil- 
ity of this system requires the satisfaction of Schwartz's 
conditions 



d^x'i" 



dx°'dx^ dx'ydx'' ' 
and by explicit calculations leads to 



(62) 



>) 

dx'^ 



-(b) 



-(a) 



{x') 



de'^'> {x) de'a' (x) 



(a) 



dx° 



dx^ 



(63) 



so that (^ becomes 

[Xa, Xb] = XaXb — XbXa = C^^b^c , (69) 

and the homogeneity is expressed as Jacobi identity 

[[Xa, Xb] , X,] + [[Xb, X,] , Xa] + [[X,, Xa] , Xfe] = 

(70) 

which, in terms of the structure constants, reads 

CibCif + C{,Cif + C{,C% = . (71) 

By dual transformation one can get the two indexes struc- 
ture constants, more convenient for calculations 



ab ~ '^abd'-^ 



(72) 



Using the properties of the tetradic basis and multiplying 
both sides of ( |63| ) by e"^^ (x) ej^^ (x) e-p {x') and some 
algebra the expression on the left hand side becomes 



where Sabc — ^"'"^ is the Levi-Civita unitary antisym- 
metric tensor (£123 — +1). The commutation rules ( |69| ) 
expressed in the new constants acquire the compact form 



M) 



[x') 



dx'^ 



He) 



{x')- 



{x') 



^abc v V r^ad \^ 
€ AfcAc — O Aw 



(73) 



deP {x') 
dx'^ 



deP jx') 
dx'^ 



Antisymmetry is implied in the definition ( fz^ ) , while the 
Jacobi identity ( [7l|) becomes 



(64) 



ebcdC^'^C''^ = 



Analogously one obtains for the right hand side an iden- 
tical expression but different only for being function of 
X. Being the transformation a; — s- x' arbitrary, both ex- 
pressions have to be equal to the same constant, giving 
the group constant of structure]^ 



(74) 



and (|6^) for the set ( |72| ) is equivalent to the vectorial 
form 



C 



ah 



ie^^^rot e 

V 



(75) 



'dP^ 
dxP 



dx- ^(-)^('>) 



= c 



ah • 



(65) 



''For the indexes of C^^;, the parentheses are unimportant. 



Any linear transformation with constant coefficients 

e(a) = ^e^") (76) 

shows the non univocal choice of the three reference vec- 
tors in the differential forms (|5j) and with respect to such 
transformations i^ab and behave like tensors. Condi- 
tion (^ is the only one to be satisfied by the structure 
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constants, considering only non equivalent combinations 
with respect to the transformation (|76|). 
The classification of non equivalent homogeneous spaces 
reduces to the determination of all non equivalent com- 
binations of the constants C"''. 
Imposing condition (Q) one gets the relations 



[-'^3, -'^i] = n2X2 + aX^ 



(77) 



where a and (77,1,712,713) are constants related to the 
structure constants. Upon reduction of all non equiv- 
alent sets under rescaling all possible uniform spaces can 
be summarized following the Bianchi classification as in 
the following table. 



Type 


a 


ni 


n2 




I 














II 





1 








VII 





1 


1 





VI 





1 


-1 





IX 





1 


1 


1 


VIII 





1 


1 


-1 


V 


1 











IV 


1 








1 


VII 


a 





1 


1 


III (a = 1) 1 

VI (a / 1) r 


a 





1 


-1 



TABLE I. Bianchi classification - Non equivalent structure 
constants 



III. PIECEWISE REPRESENTATION OF THE 
MIXMASTER 

A. Field Equations 

In a synchronous reference system, the metric for a 
homogeneous model writes 



ds^ = - Tjab (t) e^^) (x^) ej^ {x^) dx'^dx'^ 



(78) 



where the reference vectors Ca are determined through 
(|66|), once specified the structure constants. The matrix 
Vab (t) describes the temporal evolution of the tridimen- 
sional geometry, to be derived from the Einstein's field 
equations which reduce to an ordinary differential sys- 
tem, involving functions of the t only, provided the pro- 
jection of all spatial part of vectors and tensors over the 
tetradic basis chosen, using 



R. 







-Ro(a) 










= T(a)oei^^ 







(79) 



Homogeneity reflects over all scalar quantities preventing 
the presence of any spatial gradient, incompatible with 
the problem symmetry. 

The matrix rjab is the projection over the basis of the 
spatial metric 7^/3 and the role of rjab and rj""^ to raise 
and lower the indexes is clear. 

The projection of the field equations (|2l|)-(23) over the 
tetrad gives 



1 



2'*(a) 4'^(a)«^(b) 



_-L (c) 



2^ at 



dc) 



- P, 



(a) 
(&) 



(80) 
(81) 
(82) 



where 



K(a)(6) = ri{a){b) 



(83) 



and the dot is the derivative with respect to t and the 
tridimensional Ricci tensor projected 

P(a)w='?W(c)PS (84) 
in terms of the structure constants becomes 



(a)(&) 



-Cfj'^'^Cacd + C^cd^ab + ^'^cd^bcf^ ■ (85) 
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Moreover it is not necessary the explicit form of the vec- 
tor basis as functions of the coordinate for . This task is 
related to the symmetry degree of the models considered. 



B. Kasner Solution 

The simplest homogeneous cosmological model corre- 
sponds to the Bianchi type I, whose structure constants 
identically vanish, implying also the Ricci tensor compo- 
nents to be zero so that 



(86) 



Under this conditions and in empty space, equations 
(^)-(^) reduce to the system 



1 



J_d_ 



2 



(87) 



0. 



By (88) one gets the first integral 

V7«a = 2A^ = const. . 
The contraction of the a and b indexes gets 

.,a 7 2 



7 \/7 



and then 



G = const. 



(89) 



(90) 



(91) 



Without loss of generality, upon a coordinates rescaling, 
put such constant equal to one and then also 



Substituting 
constants X'i 



K = l. (92) 
|) in (E^ one finds the relation for the 



K^a — 1 



(93) 



Lowering index h in ( p9| ) gives a system of ordinary dif- 
ferential equations with respect to ^ab 



lab 



(94) 



The set of can be considered as a matrix for a linear 
transformation and by an appropriate change of coordi- 
nates {x^,x'^,x^) is reducible to a diagonal form. Given 
Pi,P2,P3 as the corresponding matrix eigenvalues, taken 
real and different, in correspondence of the normalized 
eigenvectors n'^^-', n*^^\ n^^-*, the solution of (M) is 



(3) 



(95) 



where the constant coefhcicnts for the powers of t can 
be reduced to unity, once rescaled the coordinates. If 
the tetrad vectors are parallel to the coordinate axes, say 
{x,y,z), metric reduces to 



t^P'dx^ -t^P-dy^ -t^P'dz^ 



(96) 



Constants pi,P2,P3 are three arbitrary numbers, called 
Kasner indexes, which have to satisfy the conditions 



(97) 
(98) 

Except the cases 
(0, 0, 1) and (— i, |, |), Kasner indexes have to be differ- 
ent and one of them can acquire negative value. Given 
the order 



Pi + P2 + P3 = 1 
Pi + P2 + Pa = 1 , 

as a direct consequence of (pj) and 



pi < P2 < P3 , 
the corresponding variation interval is 

4<Pi<o 

2 

< P2 < 3 
2 

3<P3<1. 



(99) 



(100) 



These numbers can also be represented in the parametric 
form 



pi (u) = 

P2 (u) = 
P3 {u) = 



1 - 


\-u 


+ 




1 + 


u 


1 - 




+ 1*2 


u 


(n 


hu) 


1 - 




+ U2 



(101) 



where the parameter u varies in the interval 1 < u < -|-oo. 
When u < 1, the parameterization can be reduced to the 
same variability interval for pi,P2,P3, holding 



Pi ( ;^ ) = Pi (u) 



P3 (u) 
P2 {u) . 



P2 - 

u 



(102) 



P3 



As functions of w, the parameters pi (it) andps (u) mono- 
tonically increase, while p2 (u) monotonically decreases. 
This metric corresponds to a flat space, even anisotropic, 
where the volume grows proportionally with increasing 
t and distances along (y, z) axes increase while along 
X decrease. The time t = is the singular point for 
the solution and such singularity is not avoidable under 
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any reference system change, while the invariant quan- 
tities of the curvature tensor diverge except the case 
Pi = P2 = 0, P3 = 1 in which the metric is reducible 
to the Galilean form once posed the parameterization 



t sinh — 
t cosh — T . 



(103) 



The metric ( pq ) is the exact solution for the Einstein's 
equations in empty space but, close the a singular point, 
say small t, it represents an approximate solution. 
A generalized solution coincides only with an approxima- 
tion of the metric, in the sense that dominant terms of 



such metric as powers of t have analogous form to ( 96 ) . 
In general, given a synchronous reference system, the 
metric can be expressed in the form ( [TtI ) , where the spa- 
tial line element dl is 

df = {a^lalp + iP'mamp -f c^UaUp) dx°'dx'^ (104) 

once posed 



c = tP" 



and the three tridimensional vectors 1, m, n redefine the 
directions along whic h the spatial lengths vary following 
the power laws (105). Such quantities are functions of 
the spatial coordinates. In view of the fact that expo- 
nents (105) are all different, the spatial metric (104) is 
always anisotropic. 

An eventual presence of matter doesn't affect the gener- 
ality of these conclusions and can be introduced in the 
metric ( |104[ )-(105) via four coordinates functions to de- 
termine the initial distribution of matter and the three 
initial velocity components. The behavior of matter in 
the vicinity of a singular point is determined by equa- 
tions of motion in a given gravitational field, following 
classical hydrodynamics. 



C. Oscillatory Regime 

Differently from the Kasner solution, in all other ho- 
mogeneous models the projection of the tridimensional 
Ricci tensor has components different from zero and the 
complexification induced prevents an analytical descrip- 
tion of the solution, except for the Bianchi type II model. 
In the following is described the asymptotic behavior of 
the homogeneous models in the vicinity of a singular 
point. 

The most general and interesting case relies in the 
Bianchi types VIII and IX models (Mixmaster) Q , while 
others can be considered as simplifications. 
The general solution is, by definition, stable. A pertur- 
bation to the system is equivalent to a change in the 



initial conditions at a given time, but the general solu- 
tion satisfies arbitrary initial conditions, then the given 
perturbation cannot affect the form of the solution. 
The structure constants sets considered are (see table) 



Bianchi VIII 
Bianchi IX : 



- ^ 1, = -1 

= = 1 = ^ 1 . 



(106) 



Let's consider the case of Bianchi IX. If the matrix rjab 
is diagonal, the components of the Ricci tensor in 
(p0|)-(|8^) vanish identically in the synchronous reference 
system and, by (|85|), the same holds for the off diagonal 
components of P{a){b)- The remaining Einstein's equa- 
tions give for the functions a{t) ,b {t) , c (t) the system 



(dbc)' 



abc 2a^ 



{b- - c^) 



{abcj 



1 



abc 2a'^b'^c? 
{aht)' _ 1 
abc 2a^6^c^ 



[a^-c^?-b^ 



2 i,2\2 4 

I — — C 



(107) 



(105) and 



a 
a 



= 



(108) 



depending on the time variable only and [)* — ■ji- Equa- 



tions (107)-(108) are exact and valid also far from the 
singularity in < = 0. Introducing the quantities 



a 



b = e' 
c = 

and the temporal variable r 

dt — abc dr 
one obtains the simplified form 



(109) 



(110) 



2arr = [b 



2 



2Prr ^ ia" ~ C^f -b^ 



-6^ - 



■{a + (3 + 7)^^ = arl3r + a^7r -I- Prlr 



(111) 



(112) 



where the index r refers to the derivative with respect to 
T. The system (111), by (112), admits the first integral 



^^(^a'^ + b^ + c^- 2a%^ - 2a\^ ~ 2b\A^ (113) 

containing only first derivatives. The Kasner regime 
(105) is a solution of ( [lllD if all right hand side expres- 
sions can be neglected. Such a situation cannot hold 
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indefinitely (for t ^ 0), because there is always an in- 
creasing term. Suppose the ca se th at, for a certain time 
interval, some expressions in ( |lll| ) are negligible (Kas- 
nerian regime): for example, the negative exponent cor- 
responds to the function a (i), so that pi — pi, and the 
perturbation to the given regime is due to terms as a*, 
meanwhil e oth er terms decrease. Keeping the dominant 
terms in (111), one reduces to the approximate 



1 4a 

arr = 

HTT ITT ^2 



(114) 



The solution for those equations describe s th e evolution 
of the metric, given the initial condition (105). Without 
loss of generality, once put 



Pi =Pi 

Pm = P2 
Pn = P3 



(115) 



so that 



a ' 
c ■ 



iP3 



(116) 



one obtains 



abc — At 
1 



A 



Int -I- const. 



(117) 



where A is a constant while initial conditions for (114) 
are 



and finally 



ar = Pi 
Pt =P2 
It =P3, 

a = Api 

P = AP2 
7 = ^P3 ■ 



(118) 



(119) 



The first of (114) has the same form as the equation 
of the unidimcnsional motion for a point particle in an 
exponential potential barrier, where a has the role of a 
coordinate. By this analogy, the initial Kasnerian regime 
corresponds to a free motion with constant velocity ar = 
pi . After the bounce on the barrier the particle follows its 
free motion with the same speed but with opposite sign 
— —pi- Using initial conditions (119), the system 



( |1 14| ) integrates as 
2 2 I pi I A 



cosh (2 I pi I At) 
= bo^ exp [2A {p2- I pi I) r] cosh (2 | pi | At) (120) 
^ co^ exp [2A (P3- I Pi I) t] cosh (2 | pi | At) 



where b p an d cq are integration constants. For the so- 
lutions (|120|) , i n the limit t ^ oo, the asymptotic form 
coincides with (llC). For t ^ ~oo, say toward the sin- 
gularity, the asymptotic forms become 



a ~ exp [Apit] 
6~exp[A(p2 + 2pi)T] 
c~exp[A(p3 + 2pi)T] 
i ~exp[A(l-H2pi)T] 



(121) 



and writing a, b, c as functions of t one gets again a Kas- 
ner behavior, a new Kasner epoch, 



b - 

c - tP'^ 
abc = A't 



(122) 



being 



Pi 



\PI 1 


1 - 


-2| 


Pi 1 


2 


Pi 


-P2 


1 


-2 


Pi\ 


P3 


- 2 


\Pi\ 


1 


- 2 


Pi\ 



(123) 



together with 



A' = (l-2|pi |)A. 



(124) 



The perturbation causes the transition from a Kasner 
epoch to another in such a way that the negative power 
of t passes from the direction 1 to m, i.e. if initially pi 
is negative, in the new solution one gets p,„' < 0. The 
initial perturbation given by terms such e^" is reduced 
substantially to zero. For the following evolution it will 
become the dominant term in another right hand side 
term and follow a completely analogous analysis. 



By the parameterization (101) and the property (102) 



the passages from a Kasner epoch to another re-express 
in a suggestive way such that if 



Pi 



- Pi {u) 

Pni = P2 (u) 
Pn = P3 (u) 



(125) 



then 



Pi 

P'm 
Pn 



= P2 (w- 1) 

= pi{u- 1) 

= P3 (u - 1) 



(126) 



and is called BKL map. 

During this transition, the function a (t) gets a maximum 
value, while b (t) a minima. After that b starts increas- 
ing, a decreases and c {t) holds its monotonic decrease. 
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During the change of Kasner epoch, the bigger of the 
t wo p ositive exponents maintains its character following 
(126). The following changes of epochs continue until 
the integer part of the initial val ue o f u becomes smaller 
then one. After that, following (102) the value of m < 1 
corresponds to another u' > 1. 

This reparameterization realizes a new exchange in the 
two positive indexes among the corresponding directions, 
say one of the oscillating functions becomes decreasing 
and vice versa. 

The time interval in which one scale factor holds mono- 
tonically decreasing while the other two oscillate is sadi 
Kasner era. The sequence of the axis exchange and the 
order in which two eraeof different time duration follows 
assumes a stochastic character. 



D. BKL Map 

Each Kasner era holds for an included number of 
epochs equivalent to the integer part of the correspond- 
ing initial value of the parameter u. Being arbitrary and 
in general irrational, the regimes' alternation continues 
indefinitely, while for an exact solution the meaning of 
the exponents Pi,P2,P3 loses the role as in the Kasner 
erae. Such internal indetermination takes out meaning to 
the study of privileged sets of values. 
The hypothesis underlying this asymptotic regime is 
based on the possibility to find a time interval during 
which it is possible to neglect the projection of the tridi- 
mensional Ricci tensor. 

The presence of matter as a perfect fluid, described by 
an opportune momentum energy tensor, doesn't affect 
the characteristics of the regime toward the singularity. 
Nevertheless the description would get a peculiar evolu- 
tion for the matter. In such a case, by the hydrodynamic 
equations, the temporal evolution of the energy density 
for each Kasner epoch becomes 



i-2(l-P3) 



(127) 



where is the bigger among the set pi,p2,P3- e mono- 
tonically increases with the temporal variable decrement 
and diverges for t — 0, as confirmation of the intrinsic 
singularity. 

In all other Bianchi models the asymptotic regime to- 
ward the singularity is overall simplified: the oscillatory 
behavio r dis appears because the right hand side mem- 
bers of (111) vanish and the corresponding Kasner epoch 
is the only one describing the system evolution toward 
the singularity. 

Oscillation amplitudes for a and f3 are growing during 
the evolution of an assigned era, independently of the era 
chosen, provided obviously the oscillating behavior. 
In parallel to that, also the Kasner epochs' duration is in- 
creasing and the matter density increases monotonically 
as 



^0 A 2k 



(128) 



being Aq the original oscillation amplitude and k the in- 
dex for the fc-th era and it is accelerating in the following 
one 



/2k' 



A 



2k^k' 



(129) 



and so on, showing how rapid is the density matter 
growth. 

The ongoing series of oscillations accumulate in the vicin- 
ity of the singular point. Between any final instant of the 
universal time t and t = Q there is an infinite number of 
oscillations. The temporal evolution is given by the nat- 
ural variable Int and not the synchronous time i, going 
to — oo toward the singularity. 

Th e study of the iterative properties of the BKL map 
( 126 ) requires another fundamental property of the pa- 
rameter u. The s-th corresponds to a succession of de- 
creasing values of the parameter u of the form u^max (star- 



mg era), u)nax 
guish 



1, U' 



II ■ 
' 1 mm' 



with notation 



u'± - < 1 

''^rriax — ^ + ^ . 



One can distin- 



(130) 



(131) 



is the maximum value of u, A;*-'''' ~ 



'^max 



anc 



x^*-* is the fractional part. The succession refers to a de- 
creasing sequence of values of u. 

The first Kasner era contains fc^**^ epochs while the sub- 
sequent one, parameterized by 



1 



(132) 



corresponds to the sequence of If the entire ar- 

rangement starts from fc^*^^ +x'^^^ , the lengths k^^^ , fc^^-* , . . . 
are the numbers involved in the expansion of x^^"^ in the 
continuous fraction 



n(0) 



1 



fc(l 



(133) 



which is finite if corresponds to the expansion of a ratio- 
nal number but in general infinite when the initial value 



The square brackets denote the integer part function. 
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is irrational. 

In any infinite sequence u constructed by ( 132 ) there are 
always arbitrarily small numbers x^"', but different from 
zero, to which correspond arbitrarily long eras k^^'^^\ 
The value assumed by all terms of such expansion is finite 
and limited and the set of all x''^^ < 1 with this property 
has null measure in the interval (0, 1). 
Depending on the strong dependence of the BKL 
map of the initial conditions the infinite sequence of 
assumes a chaotic behavior which re- 
quires an appropriate analysis. 



E. Statistical Description 

Kasner erae iteration, the specific stochasticity in the 
parameters sequence and the iterative dynamics require 
a statistical description. 

If the initial conditions don't affect the model evolution. 



expression (133) is very sensitive to the initial values of 
u and following the system along many erae it is possible 
to determine a specific probability distribution. 
Any small change on the initial value of implies a se- 
quence of numbers k completely different once explicited 
all the fraction terms, then, approaching the singularity, 
one gets a stationary probability distribution for the val- 
ues assumed by the integer part k and by the fractional 
one X referring to u. 

The random nature of the process under which the se- 
quence fc^''-' acquires asymptotically a stochastic charac- 
ter raises from the transition mechanism between differ- 



ent eras (y_32|) 



Instead of an initial value as in (13C) with s — 0, let's 
consider a distribution for cc'^"-' over the range (0, 1), de- 
scribed by the probability distribution Wq (x) for cc'"^ = 
X to be in such interval. Given Ws (x) dx the probability 
for the last term of the s-th series x^"^ = a; to be in the 
interval dx, the corresponding last term in the set has to 
be in between and necessary for the duration of 
the s-th series to be k. In this case, the probability for 
such duration is given by 



Ws [k) = / Ws-i {x)dx, 

TTk 



(134) 



and the recurrence formula relating the probability dis- 
tribution Wg+i (x) with Ws (x) is 



whose normalized solution of (|l3q) is 



w (x) 



1 



(l + x) In 2 



(137) 



Once given the probability distribution for x, the corre- 
sponding series length k is derivable as 



W{k) 



w (x) dx 



ln2 k(k 



If 



2) 



(138) 



Finally, being k and x not independent, they have to 
admit a stationary probability distribution correlated 



w Ik, x) = 



1 



{k + x) [k + X + l)\n2 



(139) 



which, given u — k + x, rewrites as a stationary distribu- 
tion for the parameter u 



w (u) 



1 



u{u+ l)ln2 



(140) 



This is the basic result for the study of the model evo- 
lution statistical properties, for the homogeneous cosmo- 
logical model considered. 

The perturbative term of the Kasner regime in the Ein- 
stein's equations is identified as the negative power law 
and the corresponding scale factor. Nevertheless, for 
u 3> 1 in any epoch, the parameterization ( |101| ) provides 
Kasner indexes in the asymptotic form 



Pi 

P2 
P3 



1 

u 
1 

u 

1 - 



(141) 



where pi,p2 have absolute values comparable and both 
close to zero. In this scheme, the perturbation is pro- 
duced equally by terms such t'^P'^ and t^^^ . Such situation 
provides equations whose solution, far from the singular- 
ity, has an oscillatory behavior. 

The obtained relations for the statistical dynamics dur- 
ing a Kasner era loose their validity when the system 
acquires a small oscillations regime, due to its strong in- 
stability. 



Ws+l (x) 



1 



xY 



1 



(135) 



The recursive definition ( |135| ) generates Ws+n {n generic 
integer) which, for increasing s, tends to a stationary 
probability distribution w (x) where the initial conditions 
are disappeared 



w (x) 



oo 



1 



-xY 



1 



IV. COVARIANT APPROACH TO THE 
MIXMASTER CHAOS 

A. The Hamiltonian Formulation 



In order to provide a Hamiltonian formulation of the 
Mixmaster dynamics, we rewrite the line element as fol- 
(136) lows 



12 



(142) 



where N{rj) denotes the lapse function, ct' are the dual 
1-forms associated with the anholonomic basis |^ and /3y 
is a traceless 3x3 symmetric matrix diag(/?ii, /322, /Jaa); 
a, iV, (3ij are functions of rj only. Parameterizing the 
matrix /3y by the usual Misner variables 



= + V3/3_ 
P22 = Ih - \/3/3- 
/333 = -2/5+ 



(143) 



the dynamics of the Mixmaster model is described by a 
canonical variational principle 



51 



with Lagrangian L 



Ldrj^O, 



(144) 



6D r 



-a + p+ + p_ 



TV 



y (a, /?+,/?_) 



(145) 



Here ()' = ^, D = det 6"+'^': 
F (a, /3_) reads 



e"^" and the potential 



(146) 



where (+) and (— ) refers respectively to Bianchi type 
VIII and IX, and the anisotropy parameters Hi (i = 
1,2,3) denote the functions El] 



H, = 

H2 = ^ 



1 13+ + \/3/3- 



3 3q: 
1 ^ /3+-y3/3- 
3 



3a 



(147) 



2/3+ 
3a 



In the limit D ^ the second three terms of the above 
potential turn out to be negligible with respect to the 



first one. Let's introduce the new (Misner-Chitre-like) 
variables 



J(r. 



(3+ = e-^'(^V^^ - 1 cos 6* 
/3_ = _ isin6'. 



(148) 



with / denoting a generic functional for m of r, 1 < ^ < 00 
and < < 27r. Then the Lagrangian (145) reads 



L = 



6D 



N 
D 



In terms of / (r) , ^ and 9 we have 



(149) 



(150) 



and since D toward the singularity, independently of 
its particular for m, in this limit / must approach infinity. 
The Lagrangian (145) leads to the conjugate momenta 



12D 



J 



df 



N \' dr' " 



_ aL _ 12D e^f 

dL 12D 2f .^2 ,^ 
89' '^^ ^ ^ 



(151) 



N 



which by a Legendre transformation make the variational 
principle (144) assume the Hamiltonian form 



6 I I + PeO' + Prr' - ) *7 = , (152) 



being 

H = 



pe {e 1) 



2AD 



P9 



24Ve^^ . (153) 



B. Reduced Variational Principle 



The dual 1-forms of the considered models are given by: 



{a = — sinh ?/ sinh 6d<l> + cosh ^pd9 
(j^ = — cosh tp sinh ^d</) + sinh ipd9 
= cosh 9d(j} + dil) 



(Bianchi IX): 



= sin 1/ sin 
(T^ = — cos ip sin 



+ cos 
+ sin ^d0 



a = cos 



+ di) 



By variating (152) with respect to TV we get the con- 
straint = 0, which solved provides 

-p, = ^ • Uadm = ^ ■ + 24^62/ (154) 
or dr 

where 



(155) 



in terms of which the variational principle (152) reduces 
to 
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S J [PiC + pee' - fUADM) dv = 0. 



(156) 



Since the equation for the temporal gauge actuaUy reads 
12D 



N{rj) 



6 T , 

'Hadm dr 



(157) 



our analysis remains fully independent of the choice of 
the time variable until the form of / and t' is not fixed. 

The variational principle (|l5q ) provides the Hamilto- 
nian equations for ^' and 9' 



Ti-ADM 

f 



{e-i)p, 

Pe 



Hadm [e - 1) 



(158) 



Furthermore can be straightforward derived the impor- 
tant relation 

diHADMl') diHADMf) , 



drj 

i.e. explicitly 

dTi-ADM 

df 



dr] 

d{nADMf') _ diHADAlf') 



=2/ 



2Ha 



df 



-24 ■ 



df 



(159) 



DM 



). (160) 



In this reduced Hamiltonian formulation, the functional 
/ (r/) plays simply the role of a parametric function of 
time and actually the anisotropy parameters Hi {i = 
1, 2, 3) are functions of the variables ^, only 



3 3e 



cos u + V3 smW 



^2 



i/3 = - + 2 

^ 3 



- - (cos 6* -VS sine 



(161) 



3C 



■ cos 9 . 



Finally, towar d the singulari ty (D — > i.e. / — > oo) by 
the expressions (146, 15C, 161), we see that 



(162) 



Since in the domain F/f , where all the Hi are simultane- 
ously greater than 0, the potential term U = e^^V can 
be modeled by the potential walls 



'^"in this study the corresponding equations for and p'g 
are not relevant. 

^^By 0() we mean terms of the same order of the inclosed 
ones. 



f/oo = Ooo {Hi (e, 9)) + Boo {H2 (e, 0)) + Qoo (fl"3 (C, 0)) (163) 

, , r +00 if a; < 
^ \ if a; > 

therefore in Th the ADM Hamiltonian becomes (asymp- 
totically) an integral of motion 



Ha 



DAI 



Ve^ + 24 • [/ ^e^ const. 



OHadm dE 



df 



df 



= 0. 



(164) 



The key point for the use of the Misner-Chitre-like vari- 
ables relies on the independence of the time variable for 
the anisotropy parameters Hi. 



C. The Jacobi Metric and the Billiard 
Representation 

Since above we have shown that asymptotically to the 
singularity (/ ^ oo, i.e. a — oo) dHADM /df = i.e . 
'Hadm ^ e = E = const., the variational principle (15£) 
reduces to 



'^J [pidi+ped9 - Edf) - 

^5 I {j)(^di + pgd9) = 0, (165) 



where we dropped the third term in the left hand side 
since it behaves as an exact differential. 

By following the standard Jacobi procedure to re- 
duce our variational principle to a geodesic one , we set 



g Pb, and by the Hamiltonian equation ( |158| ) we 



obtain the metric 



5« -^{e~ 1) 



(166) 



By these and by the fundamental constraint relation 



{e~i)pi^ 



Pe 



e-1 



= E' 



we get 



9abX X 



a'b' _ f 



E 



{{e-^)pi 



(167) 



f'E. (168) 



By the definition x"^' = = u"^, (^68|) rewrites 



(169) 
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which leads to the key relation 
drj = 



f'E 



ds . 



(170) 



Indeed expression ( |17C| ) together with p^S_' + pg0 ' = Ef 
allows us to put the variational principle ( |l65|) in the 
geodesic form 

S J f'E drj = S J Vgabu''u'>f'E ds = 

= 5 j ^jGabU°-v}> ds = (171) 

where the metric Gab = f'Egab satisfies the normaliza- 
tion condition GabU°'u'' — 1 and therefore |^ 



dr] df 



(172) 



Summarizing, in the region Th the considered dynamical 
problem reduces to a geodesic flow on a two dimensional 
Riemannian manifold described by the line element 



ds^ = E^ 



de 



1 



(e - 1) d0' 



(173) 



Now it is easy to check that the above metric has nega- 
tive curvature, since the associated curvature scalar reads 
R — therefore the point-universe moves over a neg- 

atively curved bidimensional space on which the potential 
wall ( |146D cuts the region Th- By a way completely inde- 
pendent of the temporal gauge we provided a satisfactory 
representation of the system as isomorphic to a billiard 
on a Lobachevsky plane . 



D. Invariant Lyapunov Exponent 

In order to characterize the dynamical instability of the 
billiard in terms of an invariant treatment (with respect 
to the choice of the coordinates ^, 6), let us introduce the 
following (orthonormal) tetradic basis 



1 



E 



= 0, 



1 



(174) 



Indeed the vector is nothing else than the geodesic 
field, i.e. 



'^■^ We take the positive root since we require that the curvilin- 
ear coordinate s increases monotonically with increasing value 
of /, i.e. approaching the initial cosmological singularity. 



Dv av- k I 

ds -'- 



dv 
ds 



(175) 



while the vector u;* is parallel transported along the 
geodesic, according to the equation 



ds 



dw^ 
ds 



(176) 



where by F^^ we deno ted the Christoffel symbols con- 
structed by the metric (173). Projecting the geodesic de- 
viation equation along the vector (its component along 
the geodesic field does not provide any physical infor- 
mation about the system instability), the corresponding 
connecting vector (tetradic) component Z satisfies the 
following equivalent equation 



,PZ_ 

ds^ 



Z 
E^ 



(177) 



This expression, as a projection on the tetradic basis, is 
a scalar one and therefore completely independent of the 
choice of the variables. Its general solution reads 



Z(s) 



cie' 



C2e 



ci 2 — const. 



(178) 



and the invariant Lyapunov exponent defined as [|l6| 



In 



A„ — sup lim 



2s 



in terms of the form (178) takes the value 



A,; 



E 



> 0. 



(179) 



(180) 



When the point-universe bounces against the potential 
walls, it is reflected from a geodesic to another one thus 
making each of them unstable. Though up to the limit of 
our potential wall approximation, this result shows with- 
out any ambiguity that, independently of the choice of 
the temporal gauge, the Mixmaster dynamics is isomor- 
phic to a well described chaotic system. Equivalently, in 
terms of the BKL representation, the free geodesic mo- 
tion corresponds to the evolution during a Kasner epoch 
and the bounces against the potential walls to the tran- 
sition between two of them. By itself, the positive Lya- 
punov number ( |l80| ) is not enough to ensure the system 
chaos, since its derivation remains valid for any Bianchi 
type model; the crucial point is that for the Mixmaster 
(type VIII and IX) the potential walls reduce the config- 
uration space to a compact region (Th), ensuring that 
the positivity of Au implies a real chaotic behavior (i.e. 
the geodesic motion fills the entire configuration space). 

Summarizing, our analysis shows that for any choice of 
the time variable, we are able to give the above stochas- 
tic representation of the Mixmaster model, provided the 
factorized coordinate transformation in the configuration 
space 
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(181) 



where /, a, b± denote generic functional forms of the vari- 
ables T,9,£_. 

It is worth noting that the success of our analysis, in 
showing the time gauge independence of the Mixmaster 
chaos, relies on the use of a standard ADM reduction 
of the variational principle (which reduces the system 
by one degree of freedom) and overall because, adopt- 
ing Misner-Chitre-like variables, the asymptotic poten- 
tial walls are fixed in time. The difference between our 
approach and the one presented in Q (see also for a 
critical analysis |^) consists effectively in these features, 
though in those works is even faced the problem of the 
Mixmaster chaos covariance with respect to the choice of 
generic configuration variables. 



V. STATISTICAL MECHANICS APPROACH 

A. Covariance of the Mixmaster Invariant Measure 

In order to reformulate the description of the Mixmas- 
ter stochasticity i n te rms of the Statistical Mechanics 
|l6l| , we adopt in (156) the restricted time gauge t' = 1, 
leading to the variational principle 



(182) 



In spite of this restriction, for any assigned time vari- 
able T (i.e. 77) there exists a corresponding function / (t) 
(i.e. a set of Misner-Chitre-like variables) defined by the 
(invertible) relation 



df Ti-ADM 



dr 



12D 



N (r) e- 



-2/ 



(183) 



As a consequence of the variational principle 
have again the expression (160). 

In agreement with this scheme, in the region Th where 
the potential vanishes, we have by ( |183| ) dTi.ADAi /df = 0, 
i.e. e = E = const, (by other words the ADM Hamilto- 
nian approaches an integral of motion). 

Hence the analysis to derive the invariant measure for 
the system follows the same lines presented in [^,|6| . 

Indeed we got again a representation of the Mixmaster 
dynamics in terms of a two-dimensional point-universe 
moving within closed potential walls and over a negative 
curved surface (th e Lo bachevsky plane [Q), described by 
the line element ( |173| ). Due to the bounces against the 
potential walls and to the instability of the geodesic flow 
on such a plane, the dynamics acquires a stochastic fea- 
ture. This system, admitting an "energy-like" constant 
of motion {e = E), is well-described by a microcanonical 
ensemble, whose Liouville invariant measure reads 



dg = AS{E - e) d^d9dp^dpe , A = const. (184) 

where 6 (x) denotes the Dirac function. After the natural 
positions 



Pi = 



: cos ( 



P0 



eVC^ - lsin<?!), (185) 



being < (j> < 27r, and the integration over all possible 
values of e (depending on the initial conditions, they do 
not contain any information about the system chaos) , we 
arrive to the uniform invariant measure 



du, = d£^d6d(f)-—r 



(186) 



The validity of our potential approximation is legiti- 
mated by implementing the Landau-Raichoudhury theo- 
rem 1^ near the initial singularity (placed by convention 
in T = 0, where T denotes the synchronous time, i.e. 
dT = —N{t) dr), we easily get that D vanishes monoton- 
ically (i.e. for T ^ we have dlnD/dT > 0). In terms 
of the adopted time variable r (Z) ^ /(r) 00), we 
have 



dlnD dlnDdT 



dr 



dT dr 



dlnD 
dT 



-N{t) 



(187) 



and therefore D vanishes monotonically for increasing t 
as soon as, by (183), dT/dr > 0. 

Now the key point of our analysis is tha t an y station- 
ary solution of the Liouville theorem, like ( 160 ) , remains 
valid for any choice of the time variable r; indeed in [s^ ] 
the construction of the Liouville theorem with respect to 
the variables {^,6,(j)) shows the exis tence of such proper- 
ties even for the invariant measure (191). 

We conclude remarking how, when approaching the 
singulari ty / 00 (i.e. Hadm E), the time gauge 
relation (182) simplifies to 



df ^e-2/+3Ce^ 



dT 



12 



N{T)e 



-2/. 



(188) 



in agreement with the analysis presented in ||5q| , dur- 
ing a free geodesic motion the asymptotic functions 
C (/) J ^ (/) ;</'(/)) a^re provided by the simple system 



d^ 

df 
d9 

df 



df 



sin 



cos ( 



(189) 



^ Such a theorem, within the mathematical assumptions 
founding Einsteinian dynamics, states that in a synchronous 
reference it always exists a given instant of time in cor- 
respondence to which the metric determinant vanishes 
monotonically. 
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Once getting ^ (/) as the parametric solution 
P 



sm 

/ N 1 / 1 sin^ 
/ W = 2arctanh I 



(1 



cos 



ap cos ( 



p= A/a 



a,b — const. £ 3? 



+ b 



(190) 



it reduces, for a free geodesic motion, equation (188) to 
a simple differential one for the function / (r) . 

However, the global behavior of ^ along the whole 
geodesic flow, is describ ed b y the invariant measure (191) 
and therefore relation ( |188| ) takes a stochastic character. 
If we assign one of the two functions / (r) or N (/) ana- 
lytically, the other one acquires a stochastic behavior. We 
see how the one-to-one correspondence between any lapse 
function N (77) and the associated set of Misner-Chitre- 
like variables, ensures the covariant nature with respect 
to the time-gauge of the Mixmaster universe stochastic 
behavior. 



/-^ ,dWoo , sin0 dw^ 

V? — 1 coscp—^ — h 



^ sin (j) dwoc 



0. (194) 



If now we restrict our attention to the distribution func- 
tion on the configuration space Th 



2tt 



(195) 



by ( |l94| ) we get for such restricted form the two dimen- 
sional continuity equation 

/75 7 ,dgoc , sm<j> dgoo „ 

V^cos0^ + y|f=y^ = ■ (196) 

The microcanonical solution on the whole configuration 
space 9} reads 







V {^,0}erH 



(197) 



B. Quantum Nature of the Mixmaster Chaos 

In what follows we will consider the particular case 
/(r) = T and to make clear its asymptotic nature, we 
redefine the invariant measure as follows 



dfi — Woo (Ci ^1 0) d^dOdcj) 



1 

8^ 



;d^d9d(t). 



(191) 



Summarizing, over the reduced phase space|^ {^, 6} (g) 
the distribution Woc behaves like the step-function 







'Si 



(192) 



Once performed t he tr ansformation (185) over the Hamil- 
tonian equations (158), in the asymptotic limit for which 
U — > Uoo => e = E = const., we get in Th the free 
geodesic motion [p6| 



dO sin 
d7 " 



1 cos < 



dT 



^ sin(/) 



(193) 



This dynamical system induces the stationary conti- 
nuity equation for the distribution function Woo{£,t0^4') 
describing the ensemble representation 



The main result of |^ and |5^j59[], is the proof that 
the chaos of the Bianchi IX model above outlined is an 
intrinsic feature of its dynamics and not an effect induced 
by a particular class of references: in fact the whole MCI 
formalism and its consequences can be developed in a 
framework free from the choice of a specific time gauge. 
Since this intrinsic chaos appears close enough to the Big 
Bang, we infer that it has some relations with the inde- 
terministic quantum dynamics the model performs in the 
Planckian era. This relation between quantum and deter- 
ministic chaos is searched in the sense of a semiclassical 
limit for a canonical quantization of the model [ |6ll . 

The asymptotical principle (156) describes a two di- 



mensional anholonomic Hamiltonian system, which can 
be quantized by a natural Schroedinger approach 



dip 

= Ti.ADM'ip : 



(198) 



being ip = ■i/'(T, ^,0) the wave function for the point- 



universe and, implementing Hadm (see (164)) to an op- 
eratoij^, i.e. 



Pi 



'ih- 



Pe Pe 



-ih- 



80 



(199) 



^^The only non vanishing canonical commutation relations 
are 



^'^S^ denotes the (/)-circle. 
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the equation ( |198| ) rewrites explicitly, in the asymptotic 
limit U Uoo , 



dip 



= \ e 



d 



d 



d 



Uc 

h 



1/2 



(200) 



where we took an appropriate symmetric normal order- 
ing prescription and we left Uoo to stress that the po- 
tential cannot be neglected on the entire configuration 
space {£.,0} and, being infinity out of Tjj, it requires as 
boundary condition for ip to vanish outside the potential 
walls 



(201) 



The gztanffim equation (200) is equivalent to the Wheelcr- 
DeWitt one for the same Bianchi model, once separated 
the positive and negative frequencies solutions 1 1^ , with 
the advantage that now r is a real time variable ^ Since 
the potential walls Uoo are time independent, a solution 
of this equation has the form 



(202) 



where c„ are constant coefficients and we assumed a 
discrete "energy" spectrum because the quantum point- 
universe is restricted in the finite region Th and the po- 
sition (202) in (^OOD leads to the eigenvalue problem 



1 



d_ 



1- 



d_ 



En 



Uo 



E, 



2 



(203) 



In what follows we search the semiclassical solution of 
this equation regarding the eigenvalue Eoon as a finite 
constant (i.e. we consider the potential walls as finite) 
and only at t he e nd of the procedure we will take the 
limit for Uoo ( 163 ). 

We infer that in the semiclassical limit when Ti 
and the occupation number n tends to infinity (but nh ap- 
proaches a finite value) the wave function ip„ approaches 
a function tp as 



Jim ^„ {L0)^^{t9) 



lim E'oor, — Eo 



n — >oo 



(204) 



The expression (p is taken as a semiclassical expansion up 
to the first order, i.e. 



^^This equivalence can be easily checked by taking the square 
of the operators on both sides of the equation. 



(^,0) = \MP)exp 



(205) 



where r and S are functions to be determined. 



Substituting ( |205D in (203) and separating the real from 
the complex part we get two independent equations, i.e. 



Eo 



dS_ 



1 



dS 



^2 - 1 yoe 



classical term 



(206) 



where we multiplied both sides by h and, respectively. 



d 



dS 



^«(^^W 



o(i/h) 



(207) 



In the limit h ^ the second term of ( P06D is negligible 
meanwhile the first one reduces to the Hamilton- Jacobi 
equation 



as 

'09 



Eo 



(208) 



The solution of ( |20q ) can be easily checked to be 
f f 1 ' 



'El. 



1 



di + k d9 



(209) 



where k is an inte grat ion constant. 

We observe that (|208|), through the identifications 



ds 



as 



Pe 



is reduced to the algebraic relation 

1 



Eo 



{p^d^+P0d9) , 

(210) 

(211) 



The constraint (211) is nothing more than the asymptotic 
one H^^j^j = E^= const, and can be solved setting 



dS E, 

= — , cos 



'39 



Pe = Eoo - 1 sin ( 



(212) 



where (j) G [0, 27r[ is a momentum-function related to ^ 
and 9 by the d ynam ics, (p{t) = (p {£^{t),9{t)). On the 
other hand, by (209) we get 
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1 



Pi 



to verify the compatibility of these expressions with (21^ 
we use the equations of motion (193) which provide^ 



e 



1 



(215) 



and c is a constant of integration. 

The required compatibihty comes from the identification 
k = EooC. Since 



lim 



E 
ioo 



y{^,0}erH 



(216) 



in Th, and a continuity equation ( |196| ) which has a real 
statistical meaning. 

Since referred to stationary states (pn{£,,0), the consid- 
ered semiclassical limit has to be intended in view of a 
"macroscopic" one and is not related to the temporal 
evolution of the model ifl!! . 



DEDICATION 



we see by (209) that the solution ip (^, 9) vanishes, as due 
in presence of infinit e po tential walls, outsi de T h- 
The substitution in (207) of the positions ( ^12| ) leads to 
the new equation 



ve 



T 

1 cos ffl — 



sin (j) dr 



= 0. 



(217) 



We emphasize how this equation coincides with (196), 
provided the identification r = g^; it is just this corre- 
spondence between the statistical and the semiclassical 
quantum analysis to ensure that the quantum chaos of 
the Bianchi IX model approaches its deterministic one in 
the considered limit. 

Any constant function is a solution of (|217| ) , but the nor- 
malization condition requires r — I /47r and therefore we 
finally get 



lim 

n — >oc 



io V {^,0}^Th 



(218) 

say the limit for the quantum probability distribution as 
n ^ oo and ?i — *■ associated to the wave function 

= exp |z J {p^d^ + pgde - Eoodr) | (219) 

coincides with the classical statistical distribution on the 
microcanonical ensemble. 

Though this formalism of correspondence remains valid 
for all Bianchi models, only the types VIII and IX ad- 
mit a normalizable wave function ip{^, 0), being confined 



Indeed, apart from the bounces against the potential walls. 



This paper is dedicated to the memory of Mario Im- 
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